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Abstract 

Starting from the recognition that hadrons are not produced smoothly at phase 
transition, the fluctuation of spatial patterns is investigated by finding a measure of 
the voids that exhibits scaling behavior. The Ising model is used to simulate a cross- 
over in quark-hadron phase transition. A threshold in hadron density is used to define 
a void. The dependence of the scaling exponents on that threshold is found to provide 
. useful information on some properties of the hadronization process. The complication 

', in heavy-ion collision introduces the possibility of configuration mixing, which can 

also be studied in this approach. Numerical criteria on the scaling exponents have 
been found that can be used to discriminate phase-transition processes from other 
\ hadronization processes having nothing to do with critical phenomena. 

o t 

1 Introduction 

One of the generic properties of second-order phase transition is clustering. In the case 
of a magnetic system, it means that there are regions of all sizes, in which spins point in 
the same direction, and that the probability of having a cluster of a certain size satisfies a 
scaling law. The implication of those properties for quark-hadron phase transition (PT) is 
that hadronization does not occur uniformly at the critical temperature (T c ). At any instant 
during the entire course of the hadronization process of a quark-gluon plasma system, there 
are then clusters of hadrons separated by regions of no hadrons. In the case of heavy-ion 
collisions it is the surface of the expanding cylinder that is around T c , and the clusters are to 
appear on the two-dimensional cylindrical surface. We call the non-hadronic regions between 
the clusters voids. To detect voids will be a sign of the second-order PT. In this paper we 
study the properties of the voids and discuss how to identify their properties in heavy-ion 
experiments. 

It should be understood that voids are not frozen on a plasma surface for all times. If 
one could take instantaneous pictures of the cylindrical surface at time intervals of 1 fm/c 
apart, one would see the clustering patterns to fluctuate from picture to picture. Simulation 
on a 2D lattice using the Ising model shows that at T c the clustering patterns differ from 
one configuration to another. An example of the clusters of hadrons formed can be seen 
in Fig. 2 of Ref. 0. There are regions of high hadron density, low hadron density, and 
no hadrons. A region without hadrons (i.e. a void) consists of quarks and gluons in the 
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deconfined state at a particular instant in time; they are likely to form hadrons a little later 
in the evolution process. The reason for the voids to exist is that at the critical point the 
system is torn between being in the ordered state of confinement and the disordered state of 
deconfinement. It is this tension of coexistence of the two states at T c that is responsible for 
the many interesting behaviors of critical phenomena [f2|. Thus if the plasma volume created 
in a heavy-ion collision is hot in the interior and cools to T c on the surface, a second-order 
PT, which is assumed here, would imply that the hadrons are not produced smoothly in 
time or space. Whereas hadronic clusters in the rj-(f> space may be hard to quantify, voids 
are relatively easy to define, though not trivially. Our first problem will be to characterize 
the voids. To identify them in the theoretical laboratory is, however, much simpler than to 
find them in the experimental data. The latter problem will also be addressed in this paper. 

A simple way to appreciate voids is to examine the exclusive distribution in rapidity space 
for a hadronic collision process. Since the rapidity of each produced particle can be precisely 
determined, one can calculate the rapidity gaps between each pair of neighboring particles. 
Those rapidity gaps are the ID version of the voids in 2D. An extensive consideration of the 
properties of gaps in hadronic processes is given in Ref. ||. Although it is easy to gain a 
mental picture of gaps in hadronic collisions, it is nontrivial to define a measure of voids in 
heavy-ion collisions, in which the hadronization process extends over a long period of time. 

To gain some physical insight into the fluctuation phenomenon in particle production 
at the critical point, it is helpful to consider the production of photons at the threshold of 
lasing. The physics of single-mode lasers being completely understood, their operation near 
the threshold is known to behave as in a second-order PT ||]. When the pump parameter 
is set at the threshold of lasing, the system does not continuously produce photons as a 
function of time. The photons are produced in spurts with gaps of quiescence between 
spurts. Such fluctuations have been measured and the result in terms of factorial moments 
agree with the prediction bases on the Ginzburg-Landau description of PT ||. Those 
gaps in the time series in the photo-count problem are similar to the voids in the hadron- 
count problem in heavy-ion collisions. Since lattice QCD cannot be applied efficiently to 
the simulation of spatial patterns in the PT problems, we make use of the universal features 
of the Ginzburg- Landau theory and, in particular, employ the 2D Ising model to simulate 
the hadronic clusters in the r)-(f> plane. This approach to the problem was initiated in J7|, 
where the scaling properties of cluster production were examined. More recently, an effort 
was made to find observable critical behavior in quark-hadron PT Q. Here, we use the same 
formalism to study the properties of voids and search for observable measures that can signal 
PT in heavy-ion collisions. 

2 Hadron Production in the Ising Model 

It is known through studies in lattice QCD that the nature of the phase transition depends 
on the number of flavors and the quark masses ||. For m u = = 0, the PT is first-order 
for low m s , but second-order at high m s . For nonzero m u and m^, the former region remains 
first-order (but for even smaller values of m s ), while the latter region of high m s becomes a 
cross-over. The two regions are separated by a phase boundary that is second-order and has 
possibly the Ising critical exponents ||. For realistic quark masses we are probably in the 
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region of the cross-over, which is what we shall assume in this study. 

A cross-over means that no calculable or measurable quantities and their derivatives 
undergo any discontinuity as T is varied across T c . In the 2D Ising model the situation is 
like having a small external magnetic field so that the average magnetization of the system 
varies smoothly across T c . The phase diagrams of the order parameter versus T are very 
similar for the Ising model and the QCD problem. Since lattice QCD is so much more 
difficult to study compared to the Ising model, we shall hereafter concentrate only on the 
use of the Ising model to simulate hadron production. 

Following the formalism already described in Refs. 0, 0] let us briefly summarize how 
hadron density is defined on the Ising lattice in 2D. For a lattice of size L x L with each site 
having spin cr,-, we define the spin aligned along the overall magnetization m L = J2jei 2 a j D Y 

Sj = sgn(m L )aj (1) 

where sgnimi) stands for the sign of mj. We then define the spin of a cell of size e x e at 
location i by 

Ci = s i ( 2 ) 

where Ai is the cell block of e 2 sites at i. Since q averages over all site spins in a cell, it 
should approach zero at high T for which the system is in a disordered state, and should 
approach e 2 as T — > for which the system is in an ordered state. Note that even in the 
absence of an external magnetic field, Cj approaches +e 2 at low T, never — e 2 , because of 
our definition of Sj. Thus unlike the average magnetization (m^) of the usual Ising model 
without external field, which is zero for all T > T c , the average cell spin (cj) varies smoothly 
from high to low T, similar to the behavior of (mi) in the presence of external field. 

The hadron density, being proportioned to the order parameter, can now be defined as 

p,, = Ac 2 fl( Cj ) , (3) 

where A is an unspecified factor relating the lattice spins to the number of particles in a cell. 
In any configuration q may still fluctuate from cell to cell. We identify only the positive 
Ci cells with hadron formation, and associate the hadron density with c 2 , just as the order 
parameter in the Ginzburg-Landau formalism is associated with the square of the Ising spins 
0. If (p) denotes the average density, i.e., pi averaged over all cells on the lattice and over 
all configurations, then the dependence of (p) on T is as shown in Fig. 1. It is typical of a 
cross-over, shown in Fig. 1 of Ref. O for small quark masses. The essence of that dependence 
is that (p) decreases precipitously but smoothly, as T is increased across T c , but remains 
nonzero for a range of T above T c . Evidently, we have succeeded in simulating a cross-over 
without the explicit introduction of an external field in the Ising Hamiltonian 

H = -jJ2cr l a J . (4) 

Since the dependence of (p) on T is smooth, it is a nontrivial problem to determine the 
precise value of T c . That has been done in Ref. M by examining the scaling behavior of 
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the normalized factorial moments F q . It is found that F q behaves as M Vq only at T = 2.315 
(in units of J j kg), where M is the number of bins on the lattice. Since the critical point is 
characterized by the formation of clusters of all sizes in a scale independent way, we identify 
the critical temperature at T c = 2.315. Indeed, it has been shown in Ref. |l[ that in the 
neighborhood of T c with T < T c , (p) behaves as 

(p)-(p c )oc(T c -T)\ 77 = 1.67 , (5) 

where (p c ) is (p) at T c . There exist other measures that exhibit critical behaviors, (T c — T) , 
with negative exponents. They are discussed in |]. 

3 Scaling Behavior of Voids 

We choose to work with a lattice having L = 256 and cells having e = 4. Thus the total 
number of cells on the lattice is N c = (L/e) 2 = 64 2 . We divide the lattice into bins of size S 2 
so that each bin can contain v = (S/e) 2 cells and the lattice can contain M = (L/5) 2 bins. 
The average density of hadrons in a bin is therefore 

1 v 

Pb = - ^2 Pi , (6) 

V 8=1 

where b denotes the 6th bin. Near T c , p& fluctuates from bin to bin, especially for small 5. 
We define a bin to be "empty" when 

Pb< Po , (7) 

where po is a floor level greater than zero. This criterion is chosen to eliminate the effect 
of small fluctuations on gross behavior. That is, for the purpose of defining a void, hadron 
clusters are counted only when the hadron density is above a threshold po- Bins with very 
low hadron density, i.e., where (|7|) holds, are then regarded as empty. A void is a contiguous 
collection of empty bins. Fig. 2 illustrates a pattern of voids in a configuration generated at 
T c for M = 24 2 and po = 20 (in units of A). An open square indicates an empty bin, while 
a black square contains hadrons with > p . In that configuration there are 26 voids, the 
sizes of which are 76, 35, 9, 8, 7, ■ • ■ in descending order. It should be recognized that the 
maximum hadron density that a cell can have is p max = (e 2 ) 2 = 256, so p = 20 represents a 
threshold that is less than 8% of the maximum. 

Let Vk be the size of the fcth void (in units of bins). That is, let 

V k = J20(Po~p b ) , (8) 

where (b) implies a sum over all empty bins that are connected to one another by at least 
one side; k simply labels a particular void. We can then define x k to be the fraction of bins 
on the lattice that the fcth void occupies: 

x k = Vk/M. (9) 
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For each configuration we thus have a set S = {x\, X2, ■ ■ •} of void fractions that characterizes 
the spatial pattern. 

Since the pattern fluctuates from configuration to configuration, S cannot be used to 
compare patterns in an efficient way. For a good measure to facilitate the comparison, let 
us first define the moments g q for each configuration 

-1 m 

9* = -E*l > ( 10 ) 

" l k=l 

where the sum is over all voids in the configuration, and m denotes the total number of 
voids. We then define the normalized G moments 

G q = gjgl , (11) 

which depends not only on the order q, but also on the total number of bins M. Thus by 
definition Gq = G\ = 1. This G q is defined in the same spirit as that in for rapidity gaps, 
but they are not identical because the Xk here for voids do not satisfy any sum rule. It is 
also unrelated to the G moments defined earlier [[HJ for fractal analysis. Now, G q as defined 



in Eq. (|TT|) is a number for every configuration for chosen values of q and M. With q and 
M fixed, G q fluctuates from configuration to configuration and is our quantitative measure 
of the void patterns, which in turn are the characteristic features of phase transition. 

In Fig. 3 we show the probability distribution of G q for q = 6 and M = 36 2 and po = 20 
at three different values of T. The Wolff algorithm has been used in the Monte Carlo 
simulation to reduce the correlation between configuration 0, Oj. The distribution in Fig. 3 



and other quantities to be calculated below are the results obtained using 5 x 10 3 uncorrelated 
configurations. Since the value of G q fluctuates widely from configuration to configuration 
covering a range that exceeds 3 orders of magnitude, we have plotted the distribution in 
lnG q . Since both the mean and the dispersion of lnG q shown in Fig. 3 vary significantly 
with T, and with q and M not shown in Fig. 3, it is necessary for us to search for simple 
regularities in the nature of the fluctuations of G q . 

Our first step in that search is to study the M dependence of the average of G q over all 
configurations, i.e., 

( G s> = 77 E G ? > (12) 

yV e=l 

where the superscript (e) denotes the eth event (or configuration) and Af is the total number 
of events. In Fig. 4 we show (G q ) versus M in a log-log plot for T = T c , 2 < q < 8, and 
Po = 20. We find very good linear behavior; consequently, we may write 

(G q ) oc AT* . (13) 

This scaling behavior implies that voids of all sizes occur at PT. Since the moments at 
different q are highly correlated, we expect the scaling exponent 7^ to depend on q in some 
simple way. Fig. 5 shows the dependence of 7^ on q, and we find remarkable linearity. Thus 
we may write 

lq = c + cq , (14) 
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where c = 0.8. There is no obvious reason why the (/-dependence of ^ q should be so simple. 
We should regard (|14D only as a convenient parameterization of j q that allows us to focus 
on c as a numerical description of the scaling behavior of the voids at PT. 

It is evident from Fig. 3 that studying the behavior of (G q ) extracts only a limited amount 
of information about the distribution P(G q ). The fluctuation of G q from event to event can 
be quantified by the various moments of G q 

°p« = jf£{ G PY = JdG q G>P{G q ) , (15) 

among which (G q ) corresponds only to C\ A . Instead of examining a collection of G VA for 
various values of p, we consider the derivative of C Piq at p = 1 |T!J , and define 

S « = X^U = (G q lnG q ) , (16) 



dp 



where (• • •) stands for averaging over all configurations. Despite its appearance, S q is not 
entropy, but is a measure of the fluctuations of G q , when compared with (G q ) In (G q ). 
In Fig. 6 we show the power-law behavior of S q at T c 

S q oc M a " , (17) 

where the scaling exponents a q are the slopes of the straight lines in the figure. Because of 



Eq.(|13D, {G q ) In {G q ) is not power behaved, so S q — {G q )ln{G q ) would not have a scaling 
behavior as in Eq. (|l7D . For that reason we focus on the simple properties of S q . The 
dependence of o q on q is shown in Fig. 7, where a remarkable linear behavior is found. We 
use the parameterization 

<Jq = s + sq (18) 

and find s = 0.76. 

Among the quantities that are under our control in the analysis, we have studied the 
dependences on M and q. The remaining such quantity is po, while T is beyond experimental 
control, although it can be varied in the simulation. We now study the dependence of c and 
s on po an d T, which are shown in Figs. 8 and 9. Evidently, at higher values of po the 
dependence on T are more pronounced than at po — 20. Similar behavior has been found 
for Co and sq, defined in Eqs. ([14]) and fllSj) . This result is very interesting, and provides 
a possible avenue toward learning more about the nature of PT in a realistic heavy-ion 
experiment. 

If a quark-gluon plasma is formed in a heavy-ion collision, the expanding cylinder has 
high T in the interior and low T on the surface. Our modeling has been to investigate 
the properties of hadronization on the surface. Since the PT is a smooth cross-over in the 
neighborhood of T c , and also since hydrodynamical flow can lead to local fluctuations in 
temperature and radial velocity on the surface, it is realistic to expect the hadrons to form 
in a small range of T around T c , a possibility that cannot be controlled experimentally nor 
excluded theoretically. To learn whether the hadronization takes place over a range of T, we 
suggest on the basis of Figs. 8 and 9 to use p as a device to probe the properties of the PT. 
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In an analysis of the experimental data one can use a phenomenological density threshold 
to play the role of po an d vary it in the determination of the patterns of voids. From our 
study we have learned that for a wide range of po the values of c and s are not independent of 
T. It means that if the PT occurs over a range of T so that the hadrons detected are formed 
at various T around T c even within one event, then for po in that range of nonuniform c and 
s the dependences of (G q ) and S q on M would not exhibit simple power-law behavior as in 
Figs. 4 and 6, since nonuniform values of j q and a q for any q would lead to M-dependences 
that are not simple power-laws. If, by varying p to a point that corresponds to 20 in this 
study, simple scaling behaviors can be found for (G q ) and S q with c and s having roughly 
the values 0.8 and 0.76 respectively, then we can be assured that the PT is a cross-over of 
the type that we have studied here in the Ising model. The signature for having a unique 
T c for hadronization is that the scaling behavior persists at any p and the lowest values of 
c and s are at 0.7 and 0.67, respectively, when p = 50. 

4 Configuration Mixing 

The possibility of hadronization occurring at a narrow range of temperatures discussed at 
the end of the previous section is not the only complication that may occur in a heavy- 
ion collision process. A configuration that we simulate on the Ising lattice corresponds to 
the cluster and void pattern of one instant (with uncertainty 1 fm/c) in the hadronization 
history of a plasma cylinder. The final state of a collision process registered at the detector 
is a collection of all the particles produced throughout the whole evolution process in excess 
of 10 fm/c. The clusters and voids produced at different times overlap one another in the 
integration process, resulting in a smooth spatial distribution. Thus to identify the clustering 
patterns in an experiment, it is necessary to make cuts in the phase space. Since rj and (j) are 
needed to to exhibit the 2D patterns, pt is the only remaining variable, in which cuts can 
be made. Since there is some correlation between and the evolution time, the selection 
of particles having their px lying in a very narrow interval, Apt, has the effect of selecting a 
small interval, At, in evolution time However, the correspondence is not one-to-one. In 
every Ap? interval, patches of hadrons produced at neighboring times can contribute. That 
is what we mean by configuration mixing: the experimental configuration detected in a small 
Apt interval may be a mixture of parts of configurations produced at different times, the 
latter being the pure configurations that we simulate on the Ising lattice. 

To simulate a mixed configuration, we make the following choice for definiteness. We 
divide the lattice into four quadrants. In each of the quadrants we place the corresponding 
quadrant of a new and independent configuration so that the mixed configuration consists of 
four parts of four configurations. On 5 x 10 3 such mixed configurations we then performed 
the same analysis as in the preceding section. The results are summarized in Figs. 10 and 
11 where, for comparison, the straight lines are reproduced from those in Figs. 5 and 7 for 
the pure configurations. The squares are the results for j q and a q calculated from the mixed 
configurations. Evidently, configuration mixing does not introduce any discernible deviation 
from the results of the pure configurations. The agreement being so good, we see no point 
in trying out other ways of mixing. The implication of the result is remarkable, but not 
surprising. The cluster and void patterns fluctuate so much that it does not matter whether 
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some pieces of the patterns come from different configurations, provided that the appropriate 
measure of the fluctuations is extracted. What we have extracted is the scaling behavior in 
M. The scaling exponents are then found to be independent of the configuration mixing. 

5 Conclusion 

We have shown in this paper that the study of voids can be very fruitful in finding signals 
of quark-hadron phase transition in heavy -ion collisions. The use of 2D Ising model has 
been effective in simulating a cross-over in the hadronization process. G q moments have 
been defined to quantify the dependence of the voids on the bin sizes. The scaling behavior 
that has been found provides an efficient way to use the scaling exponents 7 g and o q to 
characterize the properties of the phase transition. 

The temperature at which hadronization occurs is not under experimental control. We 
have found a way to learn whether hadronization occurs at a range of T or at a unique T. 
That is achieved by varying the density threshold po- A bin whose average density is < po 
is identified as belonging to a void. The quantity p is under the control of the analyst of 
the experimental data. If by varying p one finds that a scaling behavior can be tuned out, 
i. e. the power-law dependence on M becomes invalid for a range of po, then hadronization 
does not occur at a unique T. On the other hand, if scaling remains manifest for a range of 
Po, then there is only one temperature at which hadrons are formed. 

Even in the case where hadronization takes place in a range of T, it is possible to tune po 
to a value where strict scaling can be observed. Then 7 q and a q provide the slope parameters 
c and s that can be checked as numerical constants characteristic of the cross-over PT. Since 
there are no numerical inputs in our analysis, the values c = 0.8 and s = 0.76 are predictions 
in this study. Experimental verification of those numbers would, of course, lend significant 
support to this line of study. If the experimental numbers for c and s turn out not to have 
those values, or if there is no scaling behavior at all, one could conclude that the hadrons 
are formed without the system having gone through a phase transition. 

We have further found that the study of the scaling behavior of the void moments has the 
additional virtue of being independent of configuration mixing. That property strengthens 
the argument that a small Ap? cut in the data can provide us with a window to look 
into the hadronization process at a small At time interval where hadron clusters and voids 
are formed. The effect of randomization by the possible hadron gas in the final state on 
the scaling behavior has already been found in |l| to be minimal. Thus, we have here 
a promising procedure to investigate the properties of quark-hadron phase transition that 
should be undertaken by the heavy-ion experiments. 
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Figure Captions 

Fig. 1 Average hadron density (in units of A) versus temperature (in units of J/ks) 

Fig. 2 Spatial pattern of a configuration on the 2D lattice at T c . Open squares indicate 
voids and filled squares indicate hadrons with po = 20. 

Fig. 3 Probability distributions of lnG^ (with q — 6) at two different T. 

Fig. 4 Scaling behavior of (G q ) vs M at T c . 

Fig. 5 The dependence of r y q on q. 

Fig. 6 Scaling behavior of {S q ) vs M at T c . 

Fig. 7 The dependence of a q on q. 

Fig. 8 The dependences of the slope parameter c on p and T. 
Fig. 9 The dependences of the slope parameter s on p and T. 

Fig. 10 A comparison of j q between the mixed configurations (squares) and the pure con- 
figurations (straightline taken from Fig. 5). 

Fig. 11 A comparison of a q between the mixed configurations (squares) and the pure con- 
figurations (straightline taken from Fig. 7). 
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